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Abstract 

We give an in-depth analysis of the determination procedure of the recently 
proposed heterotic nucleon distribution amplitude which hybridizes the best 
features of the Chernyak-Ogloblin-Zhitnitsky and the Gari-Stefanis models. 
With respect to the QCD sum-rule constraints, optimized versions of these 
amplitudes are derived in terms of which a "hybridity" angle can be intro- 
duced to systematically classify all models. 
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Large-momentum transfer exclusive reactions comprise two mechanisms based on the 
factorization theorem []I],|2| : the hard-gluon exchanges that describe short-range interactions 
via perturbation theory and the long-range (nonperturbative) soft processes that ensure con- 
finement in the hadronic states. Thus, such processes depend in detail on the composition of 
the hadron wave function or more precisely, the hadron distribution amplitude: $(xj,Q 2 )- 
In a physical gauge as the light-cone gauge A + = A + A 3 = 0, this is the fundamental 
covariant amplitude describing the correlations between massless valence quarks with frac- 
tional momenta < Xi = ^- < 1 (X^ Xi = 1) and relative transverse momenta k^, p + being 
the longitudinal momentum of the incoming hadron. Then, apart from logarithmic scaling 
violations, 

<&(^,Q 2 ) = f [tfk^ixiMh (!) 



where [d 2 k ± ] = lQ7r 3 5^(Zi=i ^±)Uti 



16tt 3 



and ip(xi,\t±) is the lowest-twist Fock- 



space projection amplitude of the hadron bound state. The actual calculation of §(xi,Q 2 ) 
from QCD requires nonperturbative methods such as QCD sum rules ||, lattice gauge the- 
ory or the direct diagonalization of the light-cone hamiltonian within a discretized 
light-cone setup (see, eg., 0). 

In the context of QCD sum rules, useful theoretical constraints on the nucleon distribu- 
tion amplitude have been derived by a number of authors [|7|-§|] . They have been used to 
obtain model distribution amplitudes for the nucleon |7Hl0| and the A + (1232) isobar [|IT| --|i"3|j 



in terms of the eigenfunctions of the evolution equation [p]] represented by Appell polyno- 
mials ||14j| . However, the inevitable model purposes lead to emphasis on different facets 



of the complete picture. As a result, the Chernyak-Ogloblin-Zhitnitsky (COZ) model 
predicts R = {G^l/G^ ~ 0.5, whereas the Gari-Stefanis (GS) model |TIJ, inspired by a 
semiphenomenological analysis ||15|| , gives by construction a small R, viz., R = 0.097. 

One would like the sum rules eventually to force the choice between COZ-type and 
GS-type models on grounds of mathematical, rather than phenomenological, necessity. In 
a recent publication [|l6j (see also [|17|) we have suggested that there is a third distinct 



possibility for the nucleoli distribution amplitude, which actually hybridizes features of both 
types of models into a single mold. What is novel is the requirement that these models should 
be regarded as different aspects of a more fundamental unifying structure, we termed the 
"heterotic" solution. In the present work our interest in the model is primarily theoretic. 
We shall show that the heterotic solution is uniquely determined in the parameter space 
spanned by the expansion coefficients on the orthonormalized eigenf unctions of the evolution 
equation. With respect to a \ 2 criterion to describe the deviations from the sum rules, it 
corresponds to that local minimum which is associated with the smallest possible value of 
R, notably 0.1, compatible with the sum-rule constraints 

In ||T6| , |T7| we have given evidence that the heterotic model is the single best way to 
promote agreement between theory and the data on a variety of exclusive processes including 
elastic electron-nucleon scattering and charmonium decays into pp. 

In searching for a starting point from which to develop a credible nucleon distribution 
amplitude, we employ the ansatz described in ||18|| . Then the mixed-symmetry distribution 
amplitude $7v(xj) = V(xi) — A(xi) |^U| at fixed scale Q 2 can be reconstructed from its 



moments 

^ in2n3) = I [dx\ xTxTxT * N (xi). (2) 
Jo 

To this end, consider 

3 / a \ 

{z . p) -(n, + n 2+ n 3 ) ]J (iz ■ £) <S> N { Zi ■ p)\^ = (3) 
i=l ^ 1 ' 



where z is an auxiliary lightlike vector (z 2 = 0). To determine the moments, a short- distance 
operator-product expansion is performed at some spacelike momentum where quark-hadron 
duality is supposed to be valid. Such a computation involves correlation functions of the 
form (— q 2 = Q 2 ) 



j(»itwns,m)( gjjZ ) = i J d 4 xe iq-* < Q\T(0^ lU2n3) (O)O^ («)) 1^ > (z ' 7)77/ 
_ ( z . qj n 1 +n 2 +n 3 +m+3 j (n^n-^rri) 
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where Q^ nin2n3 ^ are appropriate [||[8| three-quark operators interpolating between the proton 
state and the vacuum. Since they contain derivatives, their matrix elements are related to 
moments of the covariant amplitudes V, A, and T: 

< n\O^ in2n3) (0)\P( P ) >= fa{ Z .p) "1+^3+1^ (mn 2 n 3 )_ (5) 

The factor (z ■ 7) 77 / serves to separate out the leading twist structure in the correlator; iV 7 
is the proton spinor, and fa denotes the "proton decay constant". 

Because $jv(xj) must be a solution of the evolution equation, it can be expressed in the 
form 



~ fry (Q 2 )\-Y-n 

® N (x u Q 2 ) = <S> as (x i )J2 B n®n(x i )(-^f^) , (6) 

n=0 a s\^ ) 

in which {$ n }^° are the eigenfunctions of the interaction kernel of the evolution equation, 
orthonormalized within a basis of Appell polynomials of degree M, and $ as (xj) = 120x1X2X3 



(see [0]). [The conventions and analytical expressions given in [18] are used.] The corre- 
sponding eigenvalues 7 n are identical with the anomalous dimensions of multiplicatively 
renormalizable baryonic operators of twist three |H| . Because the 7 n are positive fractional 
numbers increasing with n, higher terms in the expansion (6) are gradually suppressed. 

The expansion coefficients B n can be determined by means of moments inversion on 
account of the orhogonality of the eigenfunctions {$ n (xj)}: 

N f 1 

Evaluating the correlator in Eq. (4) for n% + n 2 + n 3 < 3 and m = 1, the coefficients B n for 
n = 0,1 ... ,5 (i.e., max(M) = 2) can be computed upon imposing on the moments of $at 
the sum-rule constraints estimated by COZ in ||. 

The starting point of our method is to consider the hyperspace induced by the expansion 
coefficients B n and look for solutions of the form given by Eq. (6) complying with the COZ 
sum-rule requirements. The deviations from the sum rules are parametrized by a x 2 criterion 
which accounts for the order of the moments (corresponding to the superior stability of lower 
moments relatively to the higher ones [0]). For every moment (k = 1, . . . , 18), we define 



xl = (xl ia) + Xl, {h) ) [1 - e(m k - Mr n MMr x - m„)] (8) 

with 

Xlia) = min(\Mr - m k \, \m k - M^\)N^\ (9) 

where N k = \M™ tn \ or \M™ ax \, whether m k lies on the left or on the right hand side of the 
corresponding sum- rule interval (xt Q t = SfcXfc)- 

Deviations from lower-order moments are weighted by a larger penalty factor than those 
from higher-order moments: 



100, 1< k < 3 



Xk,(b) 



10, 4 < k < 9 (10) 
1, 10 < k < 18. 



In this way we were able to filter out a COZ-type solution (designated as COZ opt ) which 
is the absolute minimum of x 2 an d thus provides the best possible agreement with the sum 
rules. Evolving that solution towards lower values of R, we generated a series of local minima 
of x 2 1 plotted in Fig. 1(a) in the (B iy R) plane. [The assignments of models to symbols 
are given in Table 1.] The dashed line is a convenient fit expressed by R = 0.436415 — 
0.005374^4 - 0.00019751. The ultimate local minimum accessible from COZ opt on that x 2 
orbit corresponds to R = 0.1 — the heterotic solution. This solution, although degenerated 
with respect to R, is distinct from the GS one. The latter, as well as its optimized versions we 
determined, constitute an isolated region (an "island") in the parameter space. Technically 
this means that they correspond to local minima of x 2 at considerably lower levels of accuracy 
so that they are separated from the COZ-Het orbit by a large x 2 barrier. Without considering 
here applications to physical processes, we only mention that the predictions extracted from 
the heterotic model are dramatically different compared to those following from the GS 



model |16| , p!7| , pT| . Also, the CZ and KS models, although in the vicinity of the COZ-Het 
orbit, are actually isolated points because they correspond to much larger x 2 values (see 
Fig. 1(b) in correspondence with Table 1). One has to exercise a certain amount of care to 



be sure that these are the only regions in the parameter space contributing to the sum rules 
at the desired level of accuracy. 

The general picture which emerges is a pattern of nucleon distribution amplitudes which 
develops into several regions of different y 2 dependence (Fig. 1(b)). The new heterotic 
amplitude and the original COZ one correspond to similar x 2 values. In contrast, the 
previous models (CZ, GS, and KS) are characterized by large deviations from the sum 
rule requirements. Optimum consistency with the sum rules is provided by the amplitude 
COZ opt , but this amplitude is not favored because it fails to predict hadronic observables 
in agreement with the existing data [21]. It is remarkable that the heterotic solution [16| 



matches the King-Sachrajda || sum-rule constraints better than the original COZ amplitude. 
This is particularly important because the KS results have been independently verified 



in |TT 



The variation with shape of the nucleon distribution amplitude with R is shown graph- 
ically, obtained by interpolating between the calculated version of the optimized COZ am- 
plitude and the heterotic one, in Fig 2. Note that although the heterotic solution was 
continuously evolved from COZ opt , its geometrical characteristics are something of a hybrid 
between COZ-type and GS-type amplitudes. [The profile of $^- et is shown in [C^,r7[].] This 



heterotic character becomes apparent by considering the corresponding covariant distribu- 
tion amplitudes V, A, and T |p0 |, shown in Fig. 3. While V^et has a symmetry pattern 



similar to that of Vqs (one main maximium), T^ et is characterized by two maxima, as Tcoz- 
The inverse heterotic combination does not belong to the COZ-Het x 2_or bit. This type of 
solution can be realized by the following expansion coefficients: Bi = 4.3025, B 2 = 1.5920, 
B 3 = 1.9675, B A = -19.6580, and B 5 = 3.3531, corresponding to $ = 24.44° (see Eq. (11)) 
and x 2 = 30.634. Remarkably, this solution yields R = 0.448, i.e., it is COZ-like (albeit in 
a different region of the parameter space). The dynamic implications of our solutions for 
various physical observables will be discussed elsewhere. 

To put these remarks on mathematical grounds, we propose to consider a classification 
scheme of the various nucleon distribution amplitudes based on the observation that the op- 
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timized versions of the COZ-type and GS-type amplitudes we derived are almost orthogonal 
to each other with respect to the weight w(xi) = $ as (a;j)/120 JHJ]. Their normalized inner 
product (COZ opt , GS opt ) yields 0.1607, which corresponds to an angle of 80.8°. Thus they 
form a quasi-orthogonal basis and can be used to classify the nucleon distribution amplitudes 
in terms of a "hybridity" angle defined by 

# = arctan ( , (11) 



where Ki(i) = (GS opt ,i) and Kzij) = (COZ opt ,i) and the index % denotes one of the 
nucleon distribution amplitudes listed in Table 1. The hybridity angle parametrizes the 
mingling of geometrical characteristics attributed to COZ-like and GS-like amplitudes 
and provides a quantitative measure for their presence in any solution conforming with 
the sum- rule constraints (Fig. 3(c)). The superimposed dashed line is a fit given by 
$/[deg] = 8.5693 + O.OI6OE4 + 0.0073^ - 0.00067B|. The mixing of different geometrical 
characteristics is particularly relevant for the heterotic solution, which generically amalga- 
mates features of both types of amplitudes. In this role, the heterotic model has the special 
virtue of simultaneously fitting the twin hopes for making reliable predictions with respect 
to the experimental data while being in agreement with the sum-rule constraints |T3],[T(|PT 



The other considered models can be tuned to fit some aspects of the data, but never all 
aspects simultaneously. 

From the above considerations we conclude that a key clue to the determination of 
nucleon distribution amplitudes on the basis of QCD sum rules is the resulting pattern of 
order effected in Fig. 1, the ordering parameters being the expansion coefficient B4 and the 
hybridity angle d. This classification scheme can valuably supplement the standard fitting 
procedure to the sum rules when enlarging the basis of eigenfunctions of the evolution kernel 
to include higher-order Appell polynomials. On the phenomenological side, fixing the value 
of the ratio R by experiment, one could use Fig. 1(a) to extract the corresponding nucleon 
distribution amplitude consistent with the sum-rule constraints. 
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TABLES 

TABLE I. Expansion coefficients B n , hybridity angle and values of the ratio R = |G^|/G^, 
in correspondence with those for \ 2 for the various nucleon distribution amplitudes discussed in 
the text. 



Model 


Bi 


B 2 


B 3 


B 4 


B 5 


tf[deg] 


R 


x 2 


Symbol 


Het. 


3.4437 


1.5710 


4.5937 


29.3125 


-0.1250 


-1.89 


.104 


33.48 


• 


COZ o P t 


3.5268 


1.4000 


2.8736 


-4.5227 


0.8002 


9.13 


.465 


4.49 


■ 


GS opt 


3.9501 


1.5273 


-4.8174 


3.4435 


8.7534 


80.87 


.095 


54.95 


▲ 


q grain 


3.9258 


1.4598 


-4.6816 


1.1898 


8.0123 


80.19 


.035 


54.11 


T 


cz 


4.3050 


1.9250 


2.2470 


-3.4650 


0.0180 


13.40 


.487 


250.07 


♦ 


coz 


3.6750 


1.4840 


2.8980 


-6.6150 


1.0260 


10.16 


.474 


24.64 


□ 


KS 


3.2550 


1.2950 


3.9690 


0.9450 


1.0260 


2.47 


.412 


116.35 


O 


GS 


4.1045 


2.0605 


-4.7173 


5.0202 


9.3014 


78.87 


.097 


270.82 


A 


Samples 





3.3125 


1.4644 


3.1438 


-1.0000 


0.8750 


7.67 


.441 


4.63 


+ 


1 


3.2651 


1.4032 


3.5466 


2.8685 


1.7954 


8.94 


.405 


5.11 


+ 


2 


3.4026 


1.4917 


3.0629 


7.3430 


0.6719 


8.75 


.385 


16.07 


+ 


3 


3.7225 


1.5030 


3.6592 


10.7265 


1.5154 


9.29 


.355 


17.78 


+ 


4 


3.8407 


1.4968 


3.2142 


14.4093 


0.8757 


10.49 


.325 


19.41 


+ 


5 


3.6544 


1.4000 


3.0993 


15.5614 


-0.1329 


6.35 


.305 


18.15 


+ 


6 


3.8607 


1.4000 


3.2375 


19.8571 


-0.1635 


6.32 


.255 


20.57 


+ 


7 


3.9783 


1.4000 


3.2706 


22.4194 


-0.4805 


5.29 


.225 


21.69 


+ 


8 


4.1547 


1.4000 


3.3756 


26.1305 


-0.5855 


5.02 


.175 


23.53 


+ 


9 


3.4044 


1.5387 


4.3094 


25.5625 


0.0625 


.01 


.153 


30.80 


+ 
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FIGURES 

FIG. 1. Classification scheme of nucleoli distribution amplitudes conforming with the sum 
rules (Table 1). (a) The ratio R = \G r l /I \/G p M as a function of the expansion coefficient B4. The 
positions of the x 2 minima are indicated, (b) Distribution of local minima of \ 2 ( on a logarithmic 
scale) plotted vs. the expansion coefficient B4. (c) Pattern of nucleon distribution amplitudes 
parametrized by the hybridity angle defined in Eq. (11), vs. the expansion coefficient B4. The 
dashed curves are the fits described in the text. 

FIG. 2. A two-dimensional example showing how the profile of the nucleon distribution ampli- 
tude $Ar changes along the COZ-Het \ 2 orbit. The amplitudes are labeled by the corresponding 
R- value. 

FIG. 3. Covariant distribution amplitudes V, A, and T associated with the optimized versions 
of the COZ and GS models in comparison with the heterotic model. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403209vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403209vl 



